Physical mechanisms governing drag reduction in turbulent Taylor-Couette flow with finite-size deformable bubbles

The phenomenon of drag reduction induced by injection of bubbles into a turbulent carrier fluid has been known for a long time; the governing control parameters and underlying physics is however not well understood. In this paper, we use three dimensional numerical simulations to uncover the effect of deformability of bubbles injected in a turbulent Taylor-Couette flow on the overall drag experienced by the system. We consider two different Reynolds numbers for the carrier flow, i.e. $Re_i=5\times 10^3$ and $Re_i=2\times 10^4$; the deformability of the bubbles is controlled through the Weber number which is varied in the range $We=0.01 - 2.0$. Our numerical simulations show that increasing the deformability of bubbles i.e., $We$ leads to an increase in drag reduction. We look at the different physical effects contributing to drag reduction and analyse their individual contributions with increasing bubble deformability. Profiles of local angular velocity flux show that in the presence of bubbles, turbulence is enhanced near the inner cylinder while attenuated in the bulk and near the outer cylinder. We connect the increase in drag reduction to the decrease in dissipation in the wake of highly deformed bubbles near the inner cylinder

Deformable ellipsoidal bubbles in Taylor-Couette flow with enhanced Euler-Lagrange tracking

In this work we present numerical simulations of $10^5$ sub-Kolmogorov deformable bubbles dispersed in Taylor-Couette flow (a wall-bounded shear system) with rotating inner cylinder and outer cylinder at rest. We study the effect of deformability of the bubbles on the overall drag induced by the carrier fluid in the two-phase system. We find that an increase in deformability of the bubbles results in enhanced drag reduction due to a more pronounced accumulation of the deformed bubbles near the driving inner wall. This preferential accumulation is induced by an increase in the resistance on the motion of the bubbles in the wall-normal direction. The increased resistance is linked to the strong deformation of the bubbles near the wall which makes them prolate (stretched along one axes) and orient along the stream-wise direction. A larger concentration of the bubbles near the driving wall implies that they are more effective in weakening the plume ejections which results in stronger drag reduction effects. These simulations which are practically impossible with fully resolved techniques are made possible by coupling a sub-grid deformation model with two-way coupled Euler-Lagrangian tracking of sub-Kolmogorov bubbles dispersed in a turbulent flow field which is solved through direct numerical simulations. The bubbles are considered to be ellipsoidal in shape and their deformation is governed by an evolution equation which depends on the local flow conditions and their surface tension

Flow induced dissolution of femtoliter surface droplet arrays

The dissolution of liquid nanodroplets is a crucial step in many applied processes, such as separation and dispersion in food industry, crystal formation of pharmaceutical products, concentrating and analysis in medical diagnosis, and drug delivery in aerosols. In this work, using both experiments and numerical simulations, we \textit{quantitatively} study the dissolution dynamics of femtoliter surface droplets in a highly ordered array under a uniform flow. Our results show that the dissolution of femoliter droplets strongly depends on their spatial positions relative to the flow direction, drop-to-drop spacing in the array, and the imposed flow rate. In some particular case, the droplet at the edge of the array can dissolve about 30% faster than the ones located near the centre. The dissolution rate of the droplet increases by 60% as the inter-droplet spacing is increased from 2.5 $\mu$m to 20 $\mu$m. Moreover, the droplets close to the front of flow commence to shrink earlier than those droplets in the center of the array. The average dissolution rate is faster for faster flow. As a result, the dissolution time $T_{i}$ decreases with the Reynolds number Re of the flow as $T_{i}\propto Re^{-3/4}$. The experimental results are in good agreement with numerical simulations where the advection-diffusion equation for the concentration field is solved and the concentration gradient on the surface of the drop is computed. The findings suggest potential approaches to manipulate nanodroplet sizes in droplet arrays simply by dissolution controlled by an external flow. The obtained droplets with varying curvatures may serve as templates for generating multifocal microlens in one array

Fluctuation-induced force in homogeneous isotropic turbulence.

Understanding force generation in nonequilibrium systems is a notable challenge in statistical physics. We uncover a fluctuation-induced force between two plates immersed in homogeneous isotropic turbulence using direct numerical simulations. The force is a nonmonotonic function of plate separation. The mechanism of force generation reveals an intriguing analogy with fluctuation-induced forces: In a fluid, energy and vorticity are localized in regions of defined length scales. When varying the distance between the plates, we exclude energy structures modifying the overall pressure on the plates. At intermediate plate distances, the intense vorticity structures (worms) are forced to interact in close vicinity between the plates. This interaction affects the pressure in the slit and the force between the plates. The combination of these two effects causes a nonmonotonic attractive force with a complex Reynolds number dependence. Our study sheds light on how length scale-dependent distributions of energy and high-intensity vortex structures determine Casimir forces

A parallel interaction potential approach coupled with the immersed boundary method for fully resolved simulations of deformable interfaces and membranes

In this paper we show and discuss the use of a versatile interaction potential approach coupled with an immersed boundary method to simulate a variety of flows involving deformable bodies. In particular, we focus on two kinds of problems, namely (i) deformation of liquid-liquid interfaces and (ii) flow in the left ventricle of the heart with either a mechanical or a natural valve. Both examples have in common the two-way interaction of the flow with a deformable interface or a membrane. The interaction potential approach (de Tullio & Pascazio, Jou. Comp. Phys., 2016; Tanaka, Wada and Nakamura, Computational Biomechanics, 2016) with minor modifications can be used to capture the deformation dynamics in both classes of problems. We show that the approach can be used to replicate the deformation dynamics of liquid-liquid interfaces through the use of ad-hoc elastic constants. The results from our simulations agree very well with previous studies on the deformation of drops in standard flow configurations such as deforming drop in a shear flow or a cross flow. We show that the same potential approach can also be used to study the flow in the left ventricle of the heart. The flow imposed into the ventricle interacts dynamically with the mitral valve (mechanical or natural) and the ventricle which are simulated using the same model. Results from these simulations are compared with ad- hoc in-house experimental measurements. Finally, a parallelisation scheme is presented, as parallelisation is unavoidable when studying large scale problems involving several thousands of simultaneously deforming bodies on hundreds of distributed memory computing processors

A mechanochemical model recapitulates distinct vertebrate gastrulation modes

During vertebrate gastrulation, an embryo transforms from a layer of epithelial cells into a multilayered gastrula. This process requires the coordinated movements of hundreds to tens of thousands of cells, depending on the organism. In the chick embryo, patterns of actomyosin cables spanning several cells drive coordinated tissue flows. Here, we derive a minimal theoretical framework that couples actomyosin activity to global tissue flows. Our model predicts the onset and development of gastrulation flows in normal and experimentally perturbed chick embryos, mimicking different gastrulation modes as an active stress instability. Varying initial conditions and a parameter associated with active cell ingression, our model recapitulates distinct vertebrate gastrulation morphologies, consistent with recently published experiments in the chick embryo. Altogether, our results show how changes in the patterning of critical cell behaviors associated with different force-generating mechanisms contribute to distinct vertebrate gastrulation modes via a self-organizing mechanochemical process.</p

A mechanochemical model recapitulates distinct vertebrate gastrulation modes

During vertebrate gastrulation, an embryo transforms from a layer of epithelial cells into a multilayered gastrula. This process requires the coordinated movements of hundreds to tens of thousands of cells, depending on the organism. In the chick embryo, patterns of actomyosin cables spanning several cells drive coordinated tissue flows. Here, we derive a minimal theoretical framework that couples actomyosin activity to global tissue flows. Our model predicts the onset and development of gastrulation flows in normal and experimentally perturbed chick embryos, mimicking different gastrulation modes as an active stress instability. Varying initial conditions and a parameter associated with active cell ingression, our model recapitulates distinct vertebrate gastrulation morphologies, consistent with recently published experiments in the chick embryo. Altogether, our results show how changes in the patterning of critical cell behaviors associated with different force-generating mechanisms contribute to distinct vertebrate gastrulation modes via a self-organizing mechanochemical process.</p